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ABSTRACT 

The problem of unstable particle decay is discussed to show how elementarity of 
a subsystem immersed in an infinitely larger environment is lost. The decay law, 
when the same kind of particles as decay product make up a thermal medium, is 
worked out in detail. The relic abundance of unstable particles does not suffer from 
the Boltzmann suppression crucial at very low temperatures, because the off shell 
contribution not considered in the Boltzmann approach, becomes dominant at low 
temperatures. The short-time behavior of the non-decay probability is also clarified, 
which is important to discuss physical relevance of the non-observation of nucleon 
decay. Two powerful methods in this respect are the operator and the path integral 
approach, both of which are reviewed. 



t Talk given at AIJIC97 on Recent Developments in Nonperturbative Quantum Field Theory , 
held at Seoul, May 26-30, 1997. To appear in the Proceedings (World Scientific). 



1 



I Introduction 



Ideal elementary processes rarely occur. They presumably never occur except 
in carefully prepared laboratory experiments. It might be useful to recall that even 
processes usually considered elementary should actually be regarded taking place in 
complex environments. This has to do with how one separates a system in question 
from a surrounding environment. For instance, the beta decay, the fundamental 
weak process n — > p + e + v e , when it occurs in nuclei, is compounded by effects 
of nuclear strong interaction with the rest of nucleons. The most spectacular of this 
kind is how nucleon decay proceeds. The extremely weak process of baryon number 
violating process at the quark-lepton level must inevitably occur in the hadronic 
environment, a nucleon. A repeated question on the nucleon decay raised in the past 
is how the extremely slow process of the lifetime of order 10 31 years or even larger is 
modified by strong interaction having the time scale of order, 10~ 23 seconds. In this 
respect one must deal with the short time limit of the decay law, and one must face 
the fact that the exponential decay law is not exact in quantum mechanics. [[I]] 

In many cases of analysis on a physical system, one focusses on a small subsystem 
put in an infinitely large closed system. Behavior of the subsystem then exhibits an 
apparent loss of elementarity due to interaction with the larger environment over 
which we have no control. Phrased somewhat differently, one often asks how a 
quantum system evolves with time under influence of an environment described by 
some mixed state. When one has a clear idea of how to separate the small system 
in question from the environment, the essential part of our problem becomes how 
elementary processes are modified in dissipative medium. Dissipation reflects the 
fact that we have ignorance of the larger environment. We believe that the problem 
should be analyzed from the first principles of quantum mechanics, at least at the 
conceptual level. In condensed matter physics the problem of this kind is generally 
known as quantum Brownian motion or quantum dissipation. 

It has become increasingly clear to us that the bulk of the past works on quantum 
dissipation relies on the simplified form of quantum friction, namely dissipation which 
is local in time. This corresponds to the exponential decay law when one examines 
the fate of some initial excitation including the unstable particle. This approximation 
is excellent in the most dominant phase of dissipation. But it fails both in early time 
and very late time behavior of the quantum system immersed in dissipative medium. 
Some formalism exists to deal with more general nonlocal dissipation, but they are 
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not very useful to the problems we would like to address. Indeed, some of our new 
results have not been recognized in the past. 

The basic idea taken by most approaches in the past |2[], and also the one I use 
here, is that existence of a continuously infinite number of environment variables 
coupled to a finite number of subsystem variables is the essential part of dissipation. 
In this view the origin of dissipation is that once something flows from the small 
system to the infinite environment, then it practically never returns to the small 
system. Dissipation, or at least something recognized as such, occurs when one does 
not and cannot make a measurement on the environment. It is well known that 
even for a pure quantum system the entropy, a measure of dissipation, is nonvanish- 
ing once one traces out a part of the quantum system. 0] Thus one postulates that 
detailed modeling of the environment and its coupling to the small system should 
be unimportant to the dissipative behavior of a small system. After integrating out 
environment variables, one should have only a few phenomenological parameters to 
describe the dissipation, ultimately obtainable from experiments. Despite this phe- 
nomenological nature one should base all auguments on rigorous quantum physics, 
and modeling is inevitable. 

The simplest, yet the most fundamental model of quantum dissipation is har- 
monic oscillator coupled to infinitely many oscillators that make up a bath in a 
mixed state. As an ideal limit of the mixed state one can also consider a pure quan- 
tum state. The nonlinearity not considered here is presumably important to restore 
the complete equilibrium of the small plus the environment system. But our interest 
is solely in the behavior of the small system, which is modified by interaction with 
a larger environment, and we are not interested in the environment part. For this 
discussion the linear approximation should be adequate. 

The problem of how a small system behaves under influence of a larger envi- 
ronment has been investigated in a variety of approaches. Two powerful methods 
to analyze this problem are the quantum Langevin equation, and the path integral 
approach. M Both approaches have merits and demerits, but when combined, they 
become very powerful. I shall first give a fundamental result from the operator ap- 
proach, which makes clear the meaning of some basic functions frequently used later. 
In the second part of this lecture I also explain the path integral approach. Our main 
results are in our two papers already in print. |5[, || 

In this lecture I shall exclusively discuss one single problem of this rather broad 
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subject. It is the decay process in medium. The problem we set up is described as 
follows. Suppose that one would like to know how the decay of unstable particles 
proceeds when they decay into two particles that also make up the bulk of the envi- 
ronment. It is not clear how the decay law known in the pure quantum case persists 
or is modified. For instance, is the well known time dependence []!]], namely the expo- 
nential law followed by the power law decay modified? An indication of substantial 
modification to this decay law is suggested by the following consideration. Some 
amount of parent particles is clearly left behind in medium, even much later than 
the decay lifetime, because in thermal environment even a heavier parent particle 
may be created by energetic daughter particles of smaller mass. If this is the case, 
what is the fraction of the parent particles left behind? We would like to answer 
these questions and elucidate the basic law by using a rigorous result of the general 
formalism. I shall first give a fairly self-contained discussion of the general formal- 
ism, because in some literature an unnecessary approximation of the local friction is 
made at the level of the general formalism, obscuring the validity of results obtained 
when the general formalism is applied to some specific problem. Our result applied 
to the unstable particle decay in thermal medium indeed casts a serious doubt on 
the familiar result based on the Boltzmann equation that uses quantities on the mass 
shell. 



II The Model 

In the present approach one models the environment by a continuously infinite 
set of harmonic oscillators of some arbitrary spectrum and couples it to the sub- 
system via a bilinear term. As will be explained below, the two particle states of 
decay product may be taken to be the environment oscillator of this kind. Let the 
subsystem variable in question be denoted by q and the environment variable by 
Q{uj). For simplicity we assume that the subsystem has one degree of freedom, but 
it should be evident to extend it to any finite number of degrees of freedom. The 
Lagrangian of our problem consists of three parts: 

L = L q [q}+L Q [Q]+L int [q,Q}. (1) 

We take for the system-environment interaction the bilinear term: 

L q = lq 2 -V(q), (2) 
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1 poo , . . 

L Q = 2l du(Q\u)-u 2 Q\u)) , (3) 

POD 

L- mt = -q / duc{u)Q{uS). (4) 



Here u c taken to be positive is the smallest of the environment frequency spectrum, 
and c(u) describes the strength distribution of environment-system interaction. In 
the present work we only consider the harmonic oscillator for the small system: 

V{q) = \uW. (5) 

The model of environment including its interaction to the subsystem is charac- 
terized by the quantity, 

r(") = -t L , (6) 

which we called the response weight. This quantity is fundamental to the rest of 
our discussion. The simplest response weight familiar in condensed matter physics 
is given in terms of three parameters, a threshold u c , an index a and a strength c, 

rH = c( W - Wc ) a . (7) 

The other hidden parameter here is the cutoff frequency, Q, above which the response 
weight vanishes, or it is nonvanishing only for uo c < uj < Q . In condensed matter 
physics the gapless case of u c = is especially popular, and the casses of a = 1, < 
1, > 1 , are called the Ohmic, sub-Ohmic, and super-Ohmic dissipation. 

Our main interest is in the unstable particle decay. In this case the form of the 
response weight is more complicated than those given above. Let us first discuss 
how one identifys the environment variable Q(u) when the same species of particles 
as the decay product make up the environment. We take a relativistic field theory 
of Yukawa type of decay interaction: 

C mt = ^x\ (8) 

where <p is the decaying parent and x the daughter particle, with \i the coupling 
of mass dimension. We shall use the capital letter M for the mass of the parent 
if and m for the mass of the daughter x- The threshold of the decay in vacuum 



is uj c = yk 2 + 4m 2 , but this is modified in thermal environment, as will be made 
clear shortly. The decaying particle, or the Fourier component of its field operator 
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q%—S d 3 x<f(x) e lk ' x , couples to the two-body operator of x% an d we identify the 
environment variable as 

Cjf (w) Qfr) = f / ^X 2 (*) e" . (9) 

The translational invariance makes each k— mode independent. It is crucial to 
realize that each momentum state of the parent particle couples to a continuously 
infinite number of two particle states of daughter particles. The approximation, 
implicit here, that two particle x states are made of independent, non- interacting 
two particle states is equivalent to that we neglect the final state interaction of x's. 
We shall omit the vector notation such that k —>■ k , and when confusion does not 
arise, we also omit the mode k altogether. 

For definiteness, we take a thermal environment of temperature T = 1/(3. The 
response weight r(u>) is then calculable using the technique of the finite temperature 
field theory. M It is the discontinuity or the imaginary part of the self-energy n(u; , k) 
of if field in thermal medium: 

r (uj) = — — f U(u - z0+) - U(u + i0 + ) ) = --Q U(u) . (10) 

Let us explain some details of this calculation of the response weight in the 
subthreshold region of \uj\ < k . The discontinuity is readily calculable from the 
imaginary-time formalism and it is given by [Q 

-9?n(w) = -£— / dE (n(E)-n(E + u)) , (11) 
lb7rfc 

uj k Am? , . 

^ = 2 ± 2V 1 -^' (12) 



where 



< E ) = ^TI ( 13 ) 



is the Planck distribution function of T = 1/(3. Since 

n{E) - n(E + u) = n(E)( 1 + n(E + u) ) - n(E + + n(E) ) , (14) 

the imaginary part ( p]) for \u\ < k is a sum of the two contributions, x + V 9 ~* X 
and its inverse process that is allowed to occur in thermal medium. Note that if can 
be off the mass shell: u 2 — k 2 ^ the if mass 2 . The factor 1 + n represents the effect 
of stimulated boson emission. 
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On the other hand, for uj > \J k 2 + 4m 2 relevant physical processes are cp <->• 
X + X- Since 



( 1 + n(E) )( 1 + n(E + w) ) - n(E)n(E +u) = l + n(E) + n(E + u) , (15) 
the imaginary part of the self-energy is given by 

-9H( W ) = ^ [ + dE(l + 2n(E)) , (16) 



/r> l l-^-- 2 (l + 2e~^). (17) 



32vr V U 1 

The last limit is valid as u — > oo. When divided by uj, the first term here T(u) = 
SII(w)/w is ~ /i 2 /(327ru;) , which is the decay rate of <p — > XX i n vacuum, including 
the effect of prolonged lifetime at w> m v . || 

The result of computation is now summarized. For uj > \/k 2 + Am? the response 
weight is || 



H 2 




r i 4m 2 2 . l-e-^+ \ 

r(w) = 32^ ^ 1 -^p + ^ ln r^^J • (18) 

For < uj < k only the second term in the bracket of Eq.(|l8D contributes. Note that 
r(u>) — > a constant (~ M/7r x decay rate in the rest frame of if) as uj — > oo. 

We note that there is a gap between k and \/k 2 + 4m 2 for which r(o;) = 0. 
Existence of the gap is important in discussing the analytic property of some basic 
functions that appear later. However, the location of the gap in thermal medium 
differs from that in vacuum. The gap ceases to exist for the massless case m — 0. 



Ill Diagonalization of the total system 

A typical question one frequently asks with regard to the system behavior is how 
a pure state of the q system evolves in time in environment described by a mixed 
state such as the thermal one of Q(u)'s. We shall first answer this by giving explicit 
time dependent operator solution for q(t) ,p(t) = q(t) that is written in terms of 
initial values, qi ,pi ,Qi(u) , P%{oj) . With this, one can clearly express correlators 
such as {q{ti)q{t2)) in terms of the probability distribution of these values in any 
initial state. Closest to our present approach is the classic work of Ullersma , with 
an important difference of the presence of a gap in the environment spectrum, which 
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yields different behaviors of the correlators. The gapless case is applied for instance 
to phonons in medium, but in the case of our interest such as the unstable particle 
decay in medium, a gap exists if the mass of daughter particle is finite. 

The model itself, a harmonic system coupled to an infinite number of harmonic 
oscillators, frequently appears in many idealized physical problems. We utilize in the 
present investigation the exact solution of the Hamiltonian eigenvalue problem to this 
system, which might also be useful in other contexts. With this exact eigen operator 
the Heisenberg time evolution becomes evident. A great advantage of this way of 
solving the present problem is that one can employ a full analogy to the scattering 
problem, especially the analyticity based on elastic unitarity. What happens is that 
an isolated spectrum of the subsystem above the two particle threshold becomes 
unstable due to the interaction with the environment, and the single particle pole 
enters into the second Riemann sheet below the cut real axis. Due to the structure 
of our model, the unitarity relation is saturated by the elastic two-body state. The 
system then becomes integrable. 

We first consider a related problem of diagonalization of the infinite dimensional 
matrix of the potential part V: 



q I, V=[ C{U)I) I (19) 



written in matrix notation. The matrix element labels u> and u>' are to be integrated 
over here. What we are doing here is not a diagonalization of the operator itself, 
which will be dealt with later. 

The potential matrix may be decomposed into the two parts as V = V + V , 
where Vo consists of the environment part alone, 

Vo= (o ^J-oo)' (20) 

> 2 |, (21) 

(22) 

I roo 

(Si\Si) = ujI 5n H 2 / duj2ujr(uj). (23) 

We may note a trivial relation: 

\Si) (<S,| = {Si\Si)Vi, V? = Vi, (24) 
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where V% is a projection operator onto a one-dimensional subspace. The fact that 
the nontrivial part V is finite dimensional is the reason this system is solvable. 

The eigenvector of diagonalized V is given with the aid of some analytic function. 
First, one defines the proper self-energy G(z) and the full propagator F(z) by 



= ..9 ■ . 2 ■ o_^.A » (26) 



1 

- z 2 + Jf+ 2nG{z) 

using an extended response weight, r(— uo) = —r(uj) for uo < —u c . The function 
F(z) has cuts along the real axis, uo > uo c and uo < —uo c . With the assumption 
of uo c > there is a gap between the two cuts. The following discontinuity relation 
holds; 

F(uo + i0 + ) - F(uo - i0 + ) = i2n r(uo)F(uo + iO + )F(uo - i0 + ) = i2ttH(uo) , (27) 

along this cut. The frequency renormalization, corresponding to the mass renormal- 
ization in field theory, is necessary when the high frequency behavior of the response 
weight is r(uo) — ► constant. The renormalized frequency is 

uo 2 R = uol + 5uo 2 = uo 2 Q -2 / duo--^ . (28) 

J u) c UO 

When this happens, one has to subtract a term in the continuous uo integral involving 
r(uo) and replace the bare uo 2 by uo\. We shall not indicate this procedure in further 
presentation, because it is fairly straightforward. 
The eigenvector of V is then given by H 



+ \uo) + c(uo)F(uo-iO + ) / duo' - \ } — \uo'). (29) 



|*(o;)) = -c(uo)F(uo |0) 

uo' 2 -uo 2 + iW 

The notation here is such that \uo) and |0) are eigenvectors of Vo: 

(Vo-uo 2 ) \uo) = 0, V o |0) = 0. (30) 

This diagonalization involves a complex phase such as 

e M^) = F ( u _ io+)/\F(u - i0+)\ . 

Acutally one can show that these are overall phases of vectors and by removing them 
one achieves diagonalization of the real symmetric matrix V by a real orthogonal 
transformation. 



9 



The matrix diagonalization precisely parallels the operator diagonalization. The 
canonical transformation from the original to the Hamiltonian eigen operator is given 
by 



QM = QM - c(u)F(u - *0+) (q - T dJ - C{J 2 \ — Q(u') ) , (31) 

P{u) = P{u) - c(u)F(co - .0+) (p-jf du' - 7T -^L p P(u/) ) , (32) 
q = - r dujc(u)F*(uj -iO + )Q(u) , (33) 

JU! C 

p = - duc{u)F*{uj-iO + )P{uj), (34) 

Jui c 

Q(co) = Q(co) + c(co) j[ y_J 2 _ i0+ J ^ ; ), (35) 

Here we find it more convenient to use 

Q(w) = e*'Q(w) 

instead of the true eigen operator by retaining those a; dependent phases. 

Assuming the canonical commutation for the original variables, one can verify 
that diagonal variables obey the correct form of the commutation relation; 

[Q(u) ,P(u')] =i5(to-Lo') 

etc. It can be proved that with the specified phases of Q(u) , P(w) the original 
variables are all hermitian, as required. 

The overlap probability of the two vectors, or equivalently the overlap between 
q and Q(u>), is 

| (0|*(w)) | 2 = c 2 H \F(u - i0 + )\ 2 = 2ujH{uj) , (37) 

H{U) = ( W 2 -wg-n(L)V + (7rr(a;)) 2 ' (38) 

with 

U(u) =V K —^- . 39 

The analytic function F(z) contains the crucial information with regard to the 
behavior of the small subsystem. The system frequency uo originally given in the 
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Lagrangian is modified by the interaction with the environment. The shifted fre- 
quency is determined by the singularity of the self-energy function F(z). With a 
condition, 

/K 2 -) > ^, (40) 
/(A) = A + KF- 1 (v / A + *0 + ), (41) 

(in this equation ~ indicates the limit from below), it can be shown that there is no 
singularity except the branch cut starting from the threshold, uj > uj c and u < — uj c , 
in the first Riemann sheet. The original pole at uo = ujq thus moves into the second 
Riemann sheet whose location is given by z obeying 

z 2 -cu 2 -2nG(z) + 2nir(z) = 0. (42) 

The imaginary part of this location gives the decay rate of any initial configuration 
of the system, as will be made evident below. 

IV Time evolution 

With the diagonal variable derived, it is easy to write the Heisenberg operator 
solution at any time t in terms of the initial operator values. For instance, 

q {t) = - dwc(u)F*(u-iO + )Q(u,t), (43) 



Q(u , t) = cosM) Qi(u) + Pi(u) . (44) 

on 

The initial values Qi , are then rewritten in terms of the original variables. 
After some straightforward calculation one finds that 



q(t) = Pi9(t) + qi g(t)- duo^o) (h*(u,t) e-^&iCwH (h.c.) ) . (45) 



oo 



p(t) = pig(t) + qi g (t)- du^fr(uj) (k*{u,t)e- i " t b i {uj) + {h.c. 



(46) 



We have introduced 

/•CO 

g(t) =2 dujH{uj) wa.(ut) (47) 
9 f°° . r(v) smjut) 

h(u,t) = f drg(r)e-^\ (49) 
Jo 

k{uj,t) = [ dr g(r)e- iUT = g(t)e~ iuJt + iujh(uj , t) . (50) 
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Furthermore, bi(uo) = ( y/u;Q(u}) + iP(uj)/y/uj )/y/2 is the annihilation operator for 
environment harmonic oscillators. 

The important time dependence is governed by the function g(t). It is useful to 
clarify the physical significance of the function g(t) in detail. For this purpose let us 
first derive a local form of the quantum Langevin equaton. Using the explicit form 
of solution, one gets after eliminating initial g« , pi dependence, 

^| + n 2 (() q + c(t) ^ = - n 2 (() /, - c(t) u - h . (si) 
nHt ) = ilfl, c(t) = ip>l, (52) 

g 9 -g 2 g 9 -g 2 

f q = r duJHu) (h*(cu,t)e- i " t b i (u;) + (h.c.)) , (53) 

f p = [ X duyffu) (k*(cu,t)e- i " t b i (u;) + (h.c.)) . (54) 

Ju) c v ' 

The quantities that determine the subsystem behavior in this equation are the time 
dependent friction (C(t)) and the time dependent frequency squared (Q 2 (t)), which 
incorporate environmental effects. Both of these are locally determined from g(t). 
Thus g(t) describes an average behavior of the system variable disregarding the 
random force from environment. On the other hand, f q and f p give the random 
force from the envrionment. 

The system variable has been determined in terms of the initial operator values 
of both the system and the environment variables. Dependence on the system initial 
values Pi , qi are given by the function g(t). Both g(t) and g(t) can be shown to obey 
integro-differential equation of the following form (y = g or g), 

^| +co 2 y + 2 jf* draj(t - r) y(r) = , (55) 

1 /'OO POO 

aj(r) = — - / duor(uo) e~ luJT = — / duo r(uo) sin(cjr) . (56) 

2 J— OO JlU c 

The two functions, g(t) and g(t), differ in their boundary conditions: g(0) = 
0, <7(0) = 1. 

The characteristic behavior of the function g(t) is that it decreases first exponen- 
tially and then finally by an inverse power of time, as can be seen in the following 
way. Using the discontinuity formula, one may rewrite the uj integration containing 

H(u) = (f(lo + i0 + ) - F(lo - i0 + ) ) I (2m) , 
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along the real axis into the F(z) integration, the complex z running both slightly 
above and below the cuts. The factor sm(ut) is replaced by ^se~ luJt in this procedure. 
A half of this complex contour can be deformed into the second sheet, and one thereby 
encounters simple poles in the second sheet. We assume for simplicity that there 
exists only a single pole in the nearby second sheet. The intregral for g(t) may then 
be expressed as the sum of the pole contribution (at z = zq with Sszq < 0) in the 
second sheet and the contribution parallel to the imaginary axis passing through 
z = uj c , both in the first (I) and in the second (II) sheet |10[|: 



g{r) = 9 (Ke~ mz ° T ) e 32 ° T 



7T 



poo 

/ dy e~ yr ( Fj (uj c + xy) - F n (u B + iy)) 
Jo 



(57) 



with K~ x = z — ttG (zo) + mr'{zo) . 

The pole contribution given by the first term describes the exponential decay 
which usually lasts very long during the most important phase of the decay period, 
while at very late times the rest of contribution gives the power law decay. In order 
to explain this late time behavior, let us take the response weight, r{uj) = c (u—u c ) a , 
for uj c < uj < Q. We assume the parameters in the range of Q ^> uo c . The late time 
behavior of g(t) is derived from the continuous part of H{u) integration and is given 
by 

2c . , cos(cj c t + £ a ) , . 

git) ^-r(« + 1) 1 c ta+l 2 ^ , (58) 

where u is the pole mass assumed to obey uj ^> the decay rate, the imaginary part 
of the pole location. Thus, the power — a — 1 of g(t) oc t~ a ~ x is related to the 
threshold behavior of the response weight. 

One can estimate the transient time from the exponential period to the power 
period by equating the two formulas of g(t) in their respective ranges, to obtain 

£ * " \ ln ( 2c r(a + l) 7 °-» ) ' (59) 

with 7 = — Qzq the decay rate. For a very small c the factor inside the logarithm 
becomes large (oc c~ 2q ~ 3 ), and by the time t* the initial population has decreased 
like 
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It may thus be claimed that the power law behavior is difficult to observe. But we 
shall later show that this may not be so in cosmology. 

We may call the approximation that neglects the non-pole contribution for g(t) 
as the resonance approximation since this approximation is equivalent to taking a 
Breit-Wigner form for H(u), 

H(u)ph-— — — , (61) 

7T (or — UJ R ) 2 + TfU} 1 

and then integrating in Eq.((P|) for g(t) in the entire range of — oo < uj < oo without 
considering the threshold effect. It may also be called the local friction approximation 
since the equation fl55p for g(t) simply reduces to 



g+uj R g + r]g = (62) 

in this case. The solution of this equation is 

g{t) « — smiuR^e-!*, (63) 

for r) <C lur. This approximation is practically very useful in many cases, but there 
are physical effects such as the final abundance of unstable particles in thermal 
medium that cannot be explained in this approximation. 

A short-time behavior of quantum dissipation is neither well described by the 
resonance approximation. We compare a precise estimate of the time dependent 
friction C(t) and the frequency Q 2 (t) to the approximate one. Assuming a smooth 
limit of the function g(t), one has 

C(t) ~ -|(^(3)(o)2_^)(o)) , (64) 
n 2 (t) ~ -^ (3) (0), (65) 

as t — >• + . On the other hand, the pole approximation, g(t) ~ 4 sin(a;t)e~ 7 * , gives 
C(t)«27, Q 2 {t) ^u 2 + 7 2 . 

More seriously, the pole approximation violates the positivity of the reduced 
density matrix: 

^- tr p 2 w 27 tr p 2 . (66) 

This relation implies that if the initial subsystem is in a pure quantum state with 
tr p 2 = 1 , it evolves into a state with tr p 2 > 1 . Because tr p = 1 , this cannot 
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be satisfied unless some diagonal element of the density matrix is negative, whose 
absolute value is larger than 1. 

A great advantage of the operator approach is that one may explicitly work out 
various correlators. For instance, 

i rti rti 

(g(*iM*2)) = - o 9(ti ~t 2 )+ dr ds g(ti - T)a R (r - s)g(t 2 - s) 
I Jo Jo 

+ 9(ti)g(t 2 ) (p 2 ) + <Kti)<7(* 2 ) (q 2 ) + ( g{h)g{t 2 ) + g^gfo) ) \ (q iPi + PlQl ) , (67) 
a R (r) = r du (2m(u) + 1) . (68) 

Jui c ZUJ 

Coincident time limits are evaluated from these, resulting in 

POO 

(q 2 (t))= / dcu(2n i (cu) + l)ricu)\h(u,t)\ 2 

Jw c 

+ g 2 {t) (p!) + g 2 (t) (qf) + g(t)g(t) (mi + , (69) 

(p 2 (t))= / du J (2n t (uj) + l)r(u)\k(uj,t)\ 2 

+ g\t) (P 2 )+ 9 2 (t)(q-) + g(t) 9 {t){p iqi + q iPi ) , (70) 
1 r°° 

- (q(t)p{t) + p(t)q(t)) = / du (2m(u;) + 1) r(u) h(u , t)k*(u , t) 

+ g(t)g(t) (p 2 ) + g(t) g (t)(q 2 ) + (g 2 (t) + g(t) 9 (PiQi + QiPi) ■ (71) 

It is an important feature of these formulas that the initial state dependence is clearly 
separated from the rest of physics as 

( 2rii(u) + 1 ) , (q 2 ) , (p?) , (p^ + q iPi ) . 



V Relic abundance of unstable particles 
in thermal medium 

We shall apply the general result thus obtained to the unstable particle decay. 
As already mentioned, in this context the system coordinate q^ refers to the Fourier 
component of the field operator of decaying particle ip. In this section we discuss the 
number operator given by 

n(t)^\( V ^ + uq 2 {t))- 1 -. (72) 
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Here the pole mass u> is chosen as the reference frequency to define the number 
operator from the Hamiltonian of the decaying field. 

We first note that the asymptotic limit of the operator is given by 

q(t) -> - [°° duJ^uj (F*{uj-iQ + ) e ~ iut b i (u) + (h.c.)) , (73) 
pit) -> % r dujujJri^j) (F*{uj-iO + )e~ iuJt b t (uj)-(h.c)) , (74) 

due to , oo) = F{u — i0 + ), which may readily be proved. Noting that 

r(co)\h(co,oo)\ 2 = H(u>) = 1 (0| ^ }) ? , (75) 
r(a;)|*(a;,oo)| 2 = co 2 H(u) = 1 1 (0|*(a;)> | 2 , (76) 

one has the asymptotic values, 

(g 2 (oc)> = r duj 1 (0| * H) |2 (2n,(c) + 1) , (77) 

Juj c zui 
roc r ; 

(p 2 (oo)} = du- I (0|*H> | 2 (2nM + 1) . (78) 

The overlap probability of the original subsystem variable with the true eigen vari- 
able of the entire system | (0|\l/(u;)) | 2 is thus fundamental to these and subsequent 
formulas. The quantity (2rij(a;) + 1) refers to the number density of the environment 
bilinear fields. We may take the value in thermal medium, 

(2n t (u) + 1) = coth(^) . (79) 

In the infinite time limit the occupation number is then 

1 f°° Sou u 2 1 

- J du coth(^-) (0 + -) H(u) - - . (80) 



n oo 



When an precise form of the overlap H used, this formula gives a 

reliable relic abundance of unstable particles. 

Let us first check that in some limit this formula gives the familiar formula for 
the abundance. When the pole term dominates, or equivalently one approximates 
H{uj) by the Breit-Wigner function, then the temperature dependent part of the 
occupation number defined by 

duo— (Q + — ) H(uo) , (81) 
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has the factor e~ u l T at low temperatures. 

In the high temperature regime the pole approximation is excellent. But this 
approximation is not good at low temperatures. Indeed, let us examine a typical 
example by taking again the form of r(uj) = c{u — u c ) a , with < a < 1 in the 
range of u c < uj < Q (Q ^> u c ) and with uj ^> Max (u c ,T) . The result is 



p « ^r(a + l)e-^T a+1 , (82) 



uj 3 



where T is the Euler's gamma function. This shows that instead of the exponential 
suppression at low temperatures what is left in medium after the decay has a power- 
law behavior of temperature dependence (oc T a+1 ). 

An implication of this behavior to the unstable particle decay, as will be made 
more explicit shortly, is that the remnant fraction in thermal medium does not suffer 
from the Boltzmann suppression factor at temperatures even much lower than the 
mass of the unstable particle. That this is possible is due to that the conventional 
approach using the approximate Boltzmann-like equation is based on S-matrix ele- 
ments computed on the mass shell, while the true quantum mechanical equation may 
contain quantities off the mass shell. As is well known, the Green's function, which 
is quantum mechanically more fundamental than the S-matrix element, does con- 
tain important contributions off the mass shell. What is called virtual intermediate 
states in elementary quantum mechanics gives rise to the off-shell contribution. Ex- 
act treatment of the problem such as ours has indeed contributions off the mass shell, 
and moreover the off-shell contribution makes up the dominant part of behaviors at 
low temperatures. 

Let us further apply these general considerations to the decay of unstable particle; 
(p — *■ x + X ■ Since we focus on the late time behavior, the initial state dependence 
disappears: in particular, whether the parent particle is or is not in thermal equilib- 
rium with the rest of medium is not important. We shall limit our discussion here 
to the decay that occurs when the parent <p becomes non-relativistic, 

w fc ~M+^l»r, (83) 

with M the <p mass. This condition is relevant in interesting cosmological problems 
of the neutron decay at the time of nucleosynthesis and GUT X boson decay at 
baryogenesis. [JTTJ 
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Since we already know the general expression for the relic abundance, what re- 
mains to be done is the mode sum over the momentum k of the unstable particle. 
First, when the pole approximation is valid, the number density is 

n « J -0- e -P(M + kV2M) = (^1)3/2 e -M/T (g4) 

This is the familiar Boltzmann suppressed formula. 

This is a bad approximation at low temperatures, T<M. To derive a precise 
formula, one has to integrate both over u and k, which is difficult to do analytically. 
As an illustration, take a constant response weight r(u) = r(oo) , hoping that the 
asymptotic region of r(oS) dominates. Then, 

1 r(oo) f°° , 2 f°° 1 
n « — r -4— - / dkk 2 duj 



2tt 2 M 3 Jo Jk e"/ T -l 

1 r(oo)T* 

= 90^3-- (85) 
Although this is not a precise calculation, it nevertheless gives a correct temperature 
dependence. 

We numerically computed all terms including the logarithmic factor in r(u) 
along with 0[m 2 ] corrections. It turns out that the total contribution is ten times 
larger than the analytic result above: in the m — > limit, 

n « ID" 2 ^ . (86) 

The main part of this large contribution comes from |a>| < k. With a dimensionless 
constant introduced by fi = gM, this gives, relative to the photon number density 



We wrote the ratio here using the numerical value Gp, the weak interaction constant 
of mass dimensions — 2 ( Gpm 2 N ~ 10~ 5 ). 

One may estimate the equal time temperature T eq at which the power contribu- 
tion becomes equal to the Boltzmann suppressed number density, to give 

^ w — — w 10~ 13 (88) 
M 30' T 3 ' 

taking as an example \x = 10 -5 M , the weak interaction strength. This number is 
in an interesting range that may affect nucleosynthesis, but we should keep in mind 
that we did not work out the relevant three body decay, n — >• p + e + v e . 
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We shall mention another application of immediate interest in cosmology; the 



heavy X boson decay for GUT baryogenesis. It has been argued |TT| that there 
exists a severe mass bound of order, 

m x > 0[a x m pl ] « 10 16 GeV , (89) 

to block the inverse process of the X boson decay so that generation of the baryon 
asymmetry proceeds with sufficient abundance of parent X particles. The usual esti- 
mate of the mass bound mentioned above is however based on the on-shell Boltzmann 
equation. More appropriate formula in this estimate is our relic number density, 

n x ^O[10- 2 }g 2 x — . (90) 
m x 

(In a more realistic estimate one should consider the X boson decay into quarks 
and leptons. But for an order of magnitude estimate difference in statistics is not 
important.) With the GUT coupling of g\j^ = 1/40, the equal temperature is 
roughly 



T eq ~—. (91) 



M 

y 

Thus, already at temperature of about half of the X mass the Boltzmann suppressed 



formula is replaced by the power formula. The kinematical condition Eq.(|89) for 
baryogenesis must be reconsidered in view of our off-shell formula. 



VI Path integral method 



The basic idea of the influence functional method [12|] is that one is interested 
in the behavior of the subsystem alone and traces out the environment variable 
altogether in the path integral formula. Furthermore, one directly deals with the 
probability instead of the amplitude. This way one can compute the reduced density 
matrix that describes a state of the small system incorporating effects of interacting 
environment. We discuss this method here, simply because this approach gives 
the non-decay probability of unstable particle which is difficult to deal with in the 
operator approach so far used. 

We define the influence functional by convoluting with the initial state of the 
environment. To do so we assume for technical reasons that initially we may take 
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an environment state uncorrected with the system. The influence functional is thus 
obtained after integrating out the environment variables: 

^[q(r),q'(r)] = J VQ(r) J VQ\r) J dQ t J dQ' t J dQ f J dQ' f 

■5(Q f - Q' f ) K ( q(r) , Q(r) ) K* ( q'(r) , Q'(r) ) p t (Q t , Q\) , (92) 

K ( q(r) , Q(t) ) = exp ( iS [Q] + ^S int [q , Q] ) , (93) 

[Q] + Sint [q,Q}= f dr ( L Q [Q] + L int [q,Q] ) . (94) 
Jo 

The influence functional is a functional of the entire path of the system q{r) and its 
conjugate path q'ij). 

Pr{Qr , Q'i) = E W n ')^n(Qi) , (0 < W n < 1) , (95) 

n 

is the initial density matrix of the environment, which can be any mixture of pure 
quantum state n with the probability w n . What deserves to be stressed is that one 
does not observe the final state of the environment, hence integration with respect 
to the final values of Qf — Q'f is performed here. 

Once the influence functional is known, one may compute the transition proba- 
bility and any physical quantities of the q— system by convoluting dynamics of the 
system under study. For instance, the transition probability is given, with introduc- 
tion of the density matrix p^ R \ by 

J dq f J dq' f r f (qf)p {R) (qf,q' f )M<l'f), (96) 
p( R ) = J Vq(r) J Vq\r) J dq t J dq[ 

MVi) HQ{r) , q'(r) ] e iS ^~ iS ^ , (97) 

where ipij's are wave functions of the initial and the final q— states, and S q [q] is the 
action of the q— system. 

The form of the influence functional is dictated by general principles such as prob- 
ability conservation and causality. Feynman and Vernon found a closed quadratic 
form consistent with these, 



Hq(r),q'(r) 



exp 



/* dr [ T ds ( i{r)a R {r - s)Z(s) + i e(r)« 7 (r - s)X(s) ) 
Jo Jo 



(98) 



with f (r) = q( T ) - q'(r) , X(r) = q(r) + q'(r) . (99) 
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Thus two real functions a^ir) are all we need to characterize the system-environment 
interaction. These are defined here in the range of r > . The fact that depends 
on the difference of time variables, r — s, is due to the assumed stationarity of the 
environment. The Feyman- Vernon formula is valid for general Lq[Q] and L q [q], not 
limited to the harmonic oscillator model if the interaction L int [q , Q] is bilinear. 

The correlation kernels appear in the influence functional as a form of the nonlocal 
interaction and they are the dissipation aj and the noise The dissipation kernel 
a j thus computed agrees with the one defined in Eq. . Let us now specialize 
to the case of the oscillator bath of temperature T = 1/(3, which is described for a 
single oscillator of frequency uj by 



Pp(Q , Q') 



\ 1/2 
UJ \ 



n coth(/3o;/2) / 

UJ 



■ cxp 



(Q 2 + Q' 2 ) cosh(/5cu) -2QQ' 



(100) 



2 sinh(^) 
The dissipation kernel is then 

Mr) = \ T duj coth(^) r(u) fT™ . (101) 

Combined together, it gives the real-time thermal Green's function: 

a(r) = a R (r) + i ai (T) = £ 4^ (ppT[Q(u k ,t) Q(u k ,0))) , (102) 

k 

oo 



a(uj} 



DC 



dra(r)e^ = iJ2 4 ( 2 \ . 6(lo 2 - uj 2 k )) . 

^ k \ u) 2 - uj? + xe eP^ - 1 ' 



k 



(103) 



As noted already, these are given in terms of the response weight r(u>), and are 
governed by the analytic function G(z). 

For the system dynamics we further assume a single harmonic oscillator of fre- 
quency ujq. In the path integral approach integration over the sum variable X(t) is 
trivial in this case, since both the local part and the nonlocal action above are linear 
in this variable: 



2 Jo 



(£(t)X(t)-ojU(t)X(t)) 

- f dr f T ds ( £(r)a fi (r - s )£( s ) + i £(r)a,(T - s)X(s) ) . (104) 
Jo Jo 
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Thus result of the path integration of the system variable X(t) gives the classical 
integro-differential equation for £(r): 



+ u%£(t) + 2 J ds^a^s-r) = 0. (105) 



f dr f T ds£{T)a R {T-s)Z{s), (106) 
Jo Jo 



dr 2 

The end result of the £ path integral then contains an integral of the form, 

rt 

10 Jo 

using the classical solution £ (r) with specified boundary conditions, £(0) = & , £ (t) = 

In the local approximation often used the dissipation kernel is taken to have the 
form of 

a I (r) = 5uj 2 5(T)+ V 5\T), (107) 

with 5u 2 representing the frequency shift and the r] term the local friction. This 
choice enables one to solve the £ equation ( |105|) by elementary means. On the other 
hand, the noise kernel is usually given by the response weight of the form, 

r(w) = ^/S), (108) 

with f(x) some cutoff function and Q a high frequency cutoff. The cutoff is needed 
to tame the high frequency integral of or(£). The simplest cutoff function f(x) = 
9(1 — x) gives an approximate form of aj(t) with the friction rj and 

5u 2 w - - r] Q . (109) 

7T 

The rest of deduction uses the Laplace transform, and we shall be brief, leaving 
technical details to our original paper. || Solution of the integro-differential equation 



(|105|) is, using g(r) defined by eq.fl47D, given as 



CM ^ 9(t-r) ( 9(t-r)g(t) \ 

with the dot denoting derivative. 

The reduced density matrix of the quantum system at any time is obtained 
from the action written in terms of the boundary values, Sd(£f ,Xf ; £j ,Xj) , by 
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convoluting with the initial density matrix of the thermal environment. This action 
is computed as 



U _ 2 V 2 _ i „ ; i 



<5 d = --^--^-wei^ + j^/^-g^^. ( m ) 

[/ = 2 /* dr [ T ds z(r) a R (r - s) z(s) , (112) 

JO JO 

2 [ dr [ dsy(T)a R (T - s)y(s) , (113) 

JO JO 

/ dr / ds (y(r)^(s) + y(s)^(r) ) a fl (r - s) , (114) 
jo jo 



V 

w 



Vir) = ^ 

Z (r) = g( t - T )- 9 ( t - r) ||, (116) 

« T) = -^-«'( 5(t - T) -^# fl )- <U7) 
For further discussion we take as the initial state a product of thermal states, 

a system of temperature To = 1//3q and an environment of temperature T = 1//3. 

We may take T = T when we apply to the decay process of excited level initially 

in thermal equilibrium. On the other hand, in the limit of T ~~ *■ it describes the 

ground state of the system harmonic oscillator. 

After a series of straightforward Gaussian integration we find the reduced density 

matrix as a function of Xf and £/, of the form, 

P {R \X f ^ f ) =2^exp[-AX 2 f -Bt 2 f + iCX f t f ], (118) 

A =h B = V h ~T^ C = h (U9) 

h = A \h(u> ,t)\ 2 ] + ^ coth(^) (g 2 + u 2 g 2 ) , (120) 

I 2 = m\ h(cu , t)k*(u ,t)] + ± coth(^) g (g +uj 2 g) , (121) 

h = X[ \k(u ,t)\ 2 } + coth(^) (g 2 +u 2 g 2 ) . (122) 

Here cu is a reference frequency taken as that of the initial system state, and equated 
here to the frequency at the pole. If one so desires, either the renormalized ur or the 
bare u Q may be taken as another choice. If we imagine a situation in which the small 
system was added to a large environment at some time, its mutual interaction being 
absent prior to the initial time, then it is appropriate to take uq as the reference 
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frequency. Since dependence on the initial state dies away quickly as time passes, 
the choice of the initial reference is not crucial for determining the behavior of states 
at late times. Both of h(u ,t) and k(u ,t) are already defined in the preceeding 
subsection. The density matrix from which any physical quantity at time t can 
be computed has explicitly been given by the discontinuity, H(u) or r(u). 

The basic quantities that appear in the reduced density matrix are related to ex- 
pectation values of the coordinate and the momentum operators at the same moment 
by 

w = u = h 2 > (123) 

(P 2 ) = 2B+^ = I 3 , (124) 

(l(gp+ pq )) = £ i = I 2 . (125) 

Thus one may write the density matrix as 



The reduced density matrix is thus characterized by expectation values of quadratic 
operators, just as in the case of pure Gaussian system without the environmental 
effect. 

It is sometimes useful to transform the density matrix in the configuration space 
to the Wigner function fw(x,p), 

/oo 
dipt R \2x^)e-^ (127) 

= ^e X p[-4.4^-<^§£>!]. (128) 

The Wigner function is expected to give the probability distribution in the phase 
space (x , p) when the semi-classical picture is valid. Expectation value of the number 
operator, namely the occupation number, in terms of the reference frequency, equated 
to the pole location oo here, is calculated most easily from the Wigner function: 
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It consists of two terms except the trivial |, the term Co/(16A) from the Gaussian 
width of the diagonal density matrix element and the rest from the kinetic term 
— j^p. This formula of course agrees with that of the previous derivation in the 
operator method. 



VII Short-time behavior of non-decay probability 

As an application of the influence functional method, I shall disucss the short- 
time behavior of the decay probability. This is an interesting problem from the 
point of the nucleon decay, as mentioned in Introduction. I shall describe some 



fundamental aspects of this problem, and leave details to our forthcoming paper. |13 
We take as the initial state of the subsystem the first excited state of harmonic 
oscillator that may be considered as a one-particle state of unstable particle: 



|i> =410), Pi(q,q') = (-) 1/2 2cgg'e-^^' 2 )/ 2 . (130) 



After some straightforward computation one obtains the reduced density matrix for 
this case: 



p R = 2 \\- ( Jo + ftX) + + ift 3 X f £ f ) exp[ - A X) - B $ + % C X f £/ 



A 

TV 

(131) 

^s77< B = ^-f>- c = k' (132) 

~' 1 - ' - N 1 2 1 t~ 1 / '2 i -2J 



J = - J[ \h(u , t) \ 2 ] , ft x = — - 2 ( g 2 + oo 2 g 2 ) , (133) 
ft = ~ ^72 ( A 2 ( 9 2 Wg 2 ) + i!(g 2 + co 2 g 2 ) - 2hh g( 9 +^g ) ) , (134) 



J; ( j2 ( V? i -.2. -.2 \ ,2, -.2 -2 .2 - ; , •., -.2 

ft - {A \K<* >t)\ 2 ]g(9 +" 2 9 ) - m h{u , t)k(u ,t)*](g 2 + u 2 g 2 



2ul{ 



(135) 



The quantities Jj are defined in Eq.( |122| ). 

The non-decay probability, or sometime called the survival probability, is defined 
as the overlap between this density matrix and that of the first excited state written 
by the final variables, q/ ,q'f- This leads to 

-5/2 

4 V v " * ' ' A' ' 4 



%^ ( (A+^-)(B + ^-) + ^- I (130) 
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(«4+f)( B + f) + x) {MB-A) + F 1 -F 2 ) 
+ ^B-A)^ 1 (B+^)+^ 2 (A+j)-^)] • (137) 

The survival probability must be derived after deviding the disconnected contribution 
that corresponds to the ground to the ground transition, 

r„ 1-1/2 



to 



/ 1 / 3 -^ + 7 + _/ 1 + -± 



4 



2u 
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For the short-time behavior this gives a leading term for the non-decay probability 
of the form, 



i-I 



duor{uo) coth(— ) \k(u ,t)\ 2 . 



(139) 



UJ J id c 2 

One is very much interested in whether the second term is proportional to time 2 or 



not. P| Naively, 



\k(uj,t)\ 2 -> g(0) 2 t 2 =t 



(140) 



as t —>■ + , but one must be careful, because the limit may not be exchangeable with 
the uj integration when the integral is conditionally convergent. 
Since 

coth(— ) = 



e P" - 1 ' 

the temperature dependent term is given by an absolutely convergent u integral, for 
which one may exchange the t — > + limit and the integral. Thus, effect dependent 
on the environment gives a term of order t 2 at short times. 

On the other hand, the environment independent term differs, depending on 
whether 

enr> 

dco r(u) 

is finite or not. If it is finite as is the case in many applications in condensed matter 
physics, the survival probability contains terms of order t 2 . But if it approaches a 
constant, r(u) — > r(oo) , then a more careful computation gives 



UJ 



(141) 
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The result thus seems to suggest that the short time limit in our boson decay 
model gives the short time limit of the exponential decay, 

e - r *«l-rt. (142) 

But, there are many details worth of further investigation which will be dealt with 



elsewhere. 13 
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